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Abstract 

Non-commutative Euclidean scalar field theory is shown to have an eigenvalue sector which 
is dominated by a well-defined eigenvalue density, and can be described by a matrix model. This 
is established using regularizations of Mg" via fuzzy spaces for the free and weakly coupled case, 
and extends naturally to the non-perturbative domain. It allows to study the renormalization of 
the effective potential using matrix model techniques, and is closely related to UV/IR mixing. 
In particular we find a phase transition for the (p^ model at strong coupling, to a phase which 
is identified with the striped or matrix phase. The method is expected to be applicable in 4 
dimensions, where a critical line is found which terminates at a non-trivial point, with nonzero 
critical coupling. This provides evidence for a non-trivial fixed-point for the 4-dimensional NC (j)'^ 
model. 
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1 Introduction 

The idea of considering quantum field theory on "quantized" or non-commutatiye spaces 
(NCFT) was put forward a long time ago [1] , and has been pursued yigorously in the past 
years; see e.g. [2, 3] for a reyiew. A yery intriguing phenomenon which was found in this 
context is the so-called UV/IR mixing [5], linking the usual UV diyergences to new sin- 
gularities in the IR. On a technical leyel, it arises because of a yery different behayior of 
planar and non-planar diagrams, which must be distinguished on NC spaces. The planar 
diagrams are essentially the same as in the commutatiye case. The non-planar diagrams 
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however lead to oscillating integrals, which are typically finite as long as the external mo- 
mentum p is non-zero, but become divergent in the limit p — > 0. This leads to serious 
obstacles to perturbative renormalization [5]. Furthermore, it appears to signal an addi- 
tional phase denoted as "striped phase" [6-8] , which arises as the minimum of the effective 
action is no longer at zero momentum. 

Because UV/IR mixing is so generic in the NC case, it is necessary to come to terms 
with it and to find suitably adapted quantization methods. The first step is clearly a 
suitable regularization of the models. This can be achieved by parametrizing the fields in 
terms of finite matrices, which is very natural on NC spaces. Several such methods are 
available by now, using e.g. using fuzzy spaces, non-commutative tori, etc. The action for 
scalar fields is then a functional of a hermitean matrix 0, where the potential TrV{(f)) looks 
like a "pure" matrix model with U (N) invariance, which is however broken by the kinetic 
term. The UV/IR mixing is expected to be recovered in the continuum limit. 

Such a regularization has been used recently to confirm numerically the non-trivial 
phase structure mentioned above in the non-commutative 0^ model [9, 10]. There has also 
been remarkable progress on the analytical side using matrix techniques: A modified 0^ 
model with an explicit IR regulator term in the action was shown to be pcrturbatively 
renormalizable [11], and certain self-dual models of NC field theory were solved exactly 
using a matrix model formulation [13]. For gauge theories, the applicability of well-known 
techniques from random matrix theory has also been shown in simple cases [14,15]. For 
the 0^ model, a similar approach using random matrix theory does not seem possible at 
first sight, lacking U{N) invariance. Nevertheless, it was conjectured in [9] that the striped 
phase should be identified with a "matrix phase'' for the fuzzy sphere, where the action 
appears to be dominated by a pure potential model in that phase. Hence a simple analytical 
approach which allows to study also scalar NCFT's with non-trivial phase structure and 
UV/IR mixing is highly desirable. In particular, it seems that the obvious parallels between 
NCFT and pure matrix models due to UV/IR mixing have not yet been fully exploited, 
apart from integrable cases [13]. 

The aim of this paper is to show that there is indeed a simple matrix model description 
which captures a certain sector of scalar NCFT, due to UV/IR mixing. This suggests a 
new approach to scalar NCFT which not only provides new insights, but also new tools to 
study the renormalization of the effective potential. The starting point is an appropriate 
parametrization for the fields: since is a hermitian matrix, it can be diagonalized as = 
U~^diag{(j)i)U where 0j are the real eigenvalues. Hence the field theory can be reformulated 
in terms of the eigenvalues (pi and the unitary matrix U. The main observation of this paper 
is now the following: the probability measure induced on the (suitably rescaled) eigenvalues 
(pi from the path integral is sharply localized, and described by an ordinary, simple matrix 
model. This means that only fields (p with a particular eigenvalue distribution characterized 
by a certain function p{s) : [~1, 1] R"^ contribute significantly to the (euclidean) path 
integral. While this is plausible using the above parametrization, it is a nontrivial statement 
which is only true in the con-commutative regime. It is established first for the free case 
(which therefore does know about non-commutativity, contrary to a common belief) and 
extends immediately to the interacting case at least on a perturbative level. This is directly 
related to UV/IR mixing, since non- planar contributions to the eigenvalue observables are 
suppressed by the oscillatory factors, while the planar contributions can be described by a 
simple matrix model without kinetic term. It is quite obvious that this will extend also to 
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the non-perturbative level, in a suitable domain. This suggests that scalar NCFT can be 
characterized by a single function p{s). 

We then work out some simple applications of this approach, which do not require 
long computations. In the weak coupling regime, this leads to a very simple method of 
computing the mass renormalization using matrix model techniques. In particular the 
standard one-loop result for the mass renormalization in the 0^ model is recovered in a 
non-standard way, and finds a natural interpretation in the matrix model. Unfortunately 
the running of the coupling does not seem to admit such a simple computation. In any 
case, we will argue that there exists a scaling limit with a non-trivial correlation length 
in the continuum limit, suggesting the existence of a renormalized 0*^ model in 2 and 4 
dimensions. 

Extending these results to the non-perturbative regime, we find a phase transition for the 
0^ model in 2 and 4 dimensions, to a phase which is tentatively identified with the striped or 
matrix phase of [6,9]. This can be compared with numerical results available for the fuzzy 
sphere in 2 dimension, with reasonable agreement which confirms the overall picture. The 
method is expected to work better in 4 dimensions, due the stronger divergences which are 
crucial for our derivation. In the 4-dimensional case, the critical line is found to terminate 
at a non-trivial point, with nonzero critical coupling Qc 0. This is expected to be a sound 
prediction, suggesting the existence of a non-trivial fixed-point for the 4-dimensional NC 
0^ model, and hence renormalizability in an appropriate sense. 

To summarize, while the dominance of planar diagrams in NCFT is very well known, 
the description of the eigenvalue sector in terms of a simple matrix model ()39|) appears to 
be new and is very practical. This provides an alternative approach to some of the results 
of [6-8] on the phase structure in NCFT, confirming the rough picture of a phase transition 
towards a phase which breaks translational invariance. However, the phase transition is 
predicted to be higher-order, as opposed to [6-8]. We find in addition a critical coupling 
(7c 7^ in 4 dimensions, and it would be very interesting to verify this numerically. 

This paper is organized as follows. Section |21 provides some background recalling the 
UV/IR mixing, and introduces the matrix regularizations used later. Section El is the 
core of this paper: after identifying the suitable observables, we show that the eigenvalue 
distribution of free fields is given by Wigner's semi-circle law, which allows to replace the 
kinetic term by the matrix model ()34j) . Interactions are included in Section 13.21 The 
corresponding reformulation of scalar NCFT using eigenvalue and angle coordinates is 
discussed in general in Section 13.31 and 13.41 including an intuitive semi-classical picture. 
This is then applied to the (p^ model in section EJ and related to some standard results for 
hermitean matrix models. These allow in Section Elto obtain the mass renormalization in a 
very simple way. The phase transition is studied in detail for 2 and 4 dimensions in section 
ini and compared with numerical results for the fuzzy sphere. We conclude in Section [7| 
with further remarks and an outlook. 

2 NC scalar fields and UV/IR mixing 

Consider scalar field theory on the non-commutative Moyal plane in even dimensions, 
with action 

S= I d'^x \]^di(l)d4 + ^mV' + |0 ^ * ^ j . (1) 
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Here is a function on M.'^, and the -k— product is the standard Moyal product of functions 
on M'^, which can be written as 

ia^b)ix) = 1^1 dSa{x+'^e.k)b{x+y)e''y , (2) 

{0'k)i = Oijkj , k-y = kiyi , 9ij = 9ji . 

This can be understood as a pull-back of the operator product of 2 operators a and b from 
a representation of the underlying (Heisenberg) algebra 

[xi, Xj] = i9ij, (3) 

using a suitable quantization map. We assume that % is non-degenerate in this paper. 

The model ((H) written down above is not well-defined as it stands, and needs regu- 
larization. The simplest way to proceed is to use a sharp UV cut-off A, which leads to 
standard computations and will be justified below. The perturbative quantization of (0) 
differs from the commutative case by the fact that planar and non-planar diagrams must 
be distinguished. The reason is that commuting 2 plane waves with wavenumbers k and 
k' produces a factor e~^^^^ , which makes non-planar loops convergent for generic external 
momenta. More explicitly, the basic one-loop planar and non-planar self-energy diagrams 
(without counting symmetry factors) are [5] 



r 



(2) 



d'^k 1 



.(2) 



/Mr, AkOp 

r^p(p) is finite as long as p ^ due to the oscillating term, but has an IR singularity 
as p ^ because the fc-integral is then divergent as usual. This is known as UV/IR 
mixing [5], and appears to be a central feature of NC field theories. It is a serious obstacle to 
perturbative renormalization, which was only overcome recently in a modified (f)'^ model [11]. 

In this paper, we shall try to turn this UV/IR mixing into a virtue, and point out that it 
is closely related to an interesting property of the scalar field in the operator formulation, 
which seems very useful and does not hold for ordinary field theories: The dynamical field 
has a well-defined eigenvalue distribution upon quantization, which is governed by a simple 
matrix model and can be studied using a saddle-point analysis. 

2.1 Matrix regularization of 

Recall that if the non-commutative algebra (jHl) is represented on a (infinite-dimensional) 
Hilbert space Ti, the integral is given by the suitably normalized trace: 

J d'^xfix) = {27iY/^Vd^ Trf = {27rey/^ Trf (5) 

where / in the rhs is the operator version of f{x) (as obtained e.g. using the Weyl quan- 
tization map). This is a manifestation of the Bohr-Sommerfeld quantization condition. 
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relating the volume of the phase space to the dimension of the Hilbert space. The last line 
holds for Sp{d)- invariant 6ij, which we assume in this paper for simplicity. 

We want to approximate this using some finite- dimensional matrix algebra. This can 
be achieved e.g. using a suitable scaling limit of the fuzzy sphere for d = 2, or more 
generally fuzzy CP" for d = 2n, see Appendix A. Indeed fuzzy CP" Mg" in a suitable 
limit, where 6ij turns out to be invariant under Sp{2n), and the propagator is the usual 
one with a sharp momentum cutoff. Another possible regularization is using so-called NC 
lattices which are products of certain ("fuzzy") NC tori [18], see also Appendix A. These are 
technically somewhat easier to handle, but lead to a modified behavior of the propagators 
for large momenta; then the most general % can be obtained. In all these regularizations, 
is a hermitian JV x Af matrix in some finite matrix algebra Mat{M ., C), and the model 

is replaced by a matrix model (where the kinetic term breaks the U{N') symmetry). In 
particular, the trace is now over the M -dimensional Hilbert-space Ti = C-^, where M is 
related to the cutoff A and 6* in a specific way ()154j) . ()118|1 depending on the regularization. 
Then V := J d'^xl = {2tt9Y/'^ Af, and we can write 

i I d'xfix) = ^ Trf. (6) 

In particular, integrals of the type J depend only on the eigenvalues of 0, and these are 
the observables we want to study. 

To make the paper most readable, the regularization using fuzzy CP"' with sharp mo- 
mentum cutoff A will be understood, while using the conventional language of Rg as much as 
possible. The results for general (non-degenerate) 9ij and somewhat modified propagators 
would be qualitatively the same. 



3 The eigenvalue distribution of the scalar field 

The basic idea is the following: Having regularized the model (P) in terms of a finite- 
dimensional hermitean matrix 0, we can diagonalize it as 

<P = U-\<Pi)U (7) 

where f/ is a unitary A/" x A/" matrix, and (0j) = diag{(j)i, ...,(j)j\f) is diagonal with real 
eigenvalues. The integration measure in the path integral can now be written as J T>(f)e~^ = 
J (i0jA^(0j) / dUe~^, where A^(0j) = Y[i<ji4>i~ 4>j)'^ is the Vandermonde-determinant and 
dU the Haar measure for SU{M'). We are interested in the probability measure or effective 
action for these eigenvalues, induced by this path integral. For this purpose, consider e.g. 
the expectation values 

( j d^xf'^ix)) = ^ j P0exp(-S)( j rf'^x02"(x)). (8) 

They are strongly divergent normally but make sense in the regularized (matrix) case; in 
particular, we do not want to replace the 0^"(a;) by some renormalized objects but simply 
keep track of their dependence on the cutoff. Since they depend only on the eigenvalues of 
the field in the matrix representation, we can determine the effective eigenvalue distri- 
bution by studying such observables. This turns out to be non-trivial already in the free 
case: 
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3.1 The free case 

We compute the observables (jH)) for g = with a sharp UV - cutoff A using Wicks theorem. 
This involves in general planar and non-planar diagrams. The simplest case is^ 



=: c{m,K)V k'^-^. (9) 

This formal result will be fully justified in Appendix A using regularizations of M^^^" in 
terms of fuzzy CP" p53|) or fuzzy tori. Here V denotes the regularized volume of M^, 
and c = c(m, A) is of order 1 for (i > 3, and c = O(lnA) for d = 2. More precisely, in 4 
dimensions we have 

and in 2 dimensions 

r Y /-A/m y Y A2 

{ d'<p\x)) = - du^— = -Hl + —) (11) 

(the subleading behavior is modified for the regularization using NC tori, see Appendix A). 
Therefore 

,(,,„,A, = |«t(l-#'"(l + fi)). ^ = 4 

Next, consider 

( ! d^X<i){xf) = ( / d''x<P{xf) Planar + ( / rf"a;</.(x)^);Von 



Planar 



which is obtained by summing over all complete contractions of a vertex with 4 legs. There 
are 2 planar and one non-planar such contractions, the latter being given by joining the 
external legs of the non-planar self-energy diagram Q. The planar contribution is simply 

( 1 d'x4>{xf)piar.ar = 2V A"^'-'\ (13) 

(2) 

since Tp (p) is independent of p. On the other hand, the non-planar contribution is sublead- 
ing (this will be discussed in detail below), since T)^p{p) is finite except for the singularity 
at p — i> 0. This is clearly related to UV/IR mixing, even though we are considering only 
the free case up to now. We therefore expect 

a f d'^x(f)(x)^) . 1 . 

' -2 +0{-—). (14) 



(^/d%0(x)2)2 'cA^-2^- 



^we basically assume that 0ij is non-degenerate in this paper, which implies that d is even. However, 
some of the results extend to the degenerate case, which will be pursued elsewhere 
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More generally, consider 



/ d'^X(t){xy^) _ / d'^X(t){xy^) Planar / d-^XcPixf^) Non-Planar 



2\n 



The first term is of order one, and simply counts the number Npianar{'2n) of planar con- 
tractions of a vertex with 2n legs. The non-planar contributions always involve oscillatory 
integrals, and do not contribute to the above ratio in the large A limit. Let us discuss them 
in more detail: Assume that the cut-off is much larger than the NC scale, 

A^e > 1, (15) 

which will be understood from now on. By rescaling the momenta k' — k/A, the above 
ratios for d > 3 have the form 

^ _ {vJd''x(l>{xf^^)Non~Planar ^ f d% JfX^iA^EKek'. ^ Q aS A 



(16) 

which vanish for large A due to the rapidly oscillating exponential. The integrals are over 
the unit ball resp. hypercube. This is established more carefully in Appendix B, where we 
show that Rnp = 0{1/A) (at least) for d > 3. In 2 dimensions, these considerations are 
more delicate as the above ratio Rnp vanishes only logarithmically, and divergences may 
also arise in the IR. One must then assume furthermore that e.g. — is fixed while A ^ oo, 
where 

ml \ (17) 

is the non-commutative mass scale. Due to this weaker logarithmic behavior, the fiuctua- 

tions arc expected to be larger for d = 2 than for d = A. We refrain here from more precise 
estimates which should be made in the context of the regularized models, see Appendix B. 
With all these assumptions, we conclude that 



/ d''x(l){xf^) _{^J d^X(t>{xy^)pia,,ar 

J ddxcl>{x)^)- J d'^xcj>{x] 



2\n 



Npianar{2n) (18) 



in the large A limit. Notice that this result is very different from the conventional field 
theory: the above calculation with a naive cutoff would have the same contributions from 
planar and non-planar diagrams. Then the total number of contractions of e.g. (0^"') is of 
order 2"n! 3> -/Vpzanar(2n), which would invalidate the conclusions below. 

Next, consider expectation values of products ((J d'^x<p{x)'^"-^)...{J d!^x(p{xf-"'')) . We 
claim that this factorizes in the large A limit, 

((^ / d'^x0(x)2^'i)...(^ / d'^x(l){xY''^)) J d'^x(t){xf"^) / d^a;0(x)2"fc) 

/ d'^x(l){xY)'^K..{^ J d'^x(l){xfY'k ^ JfJ¥^^ix)^ "' {^J ddx(l){xy)^'' 

= N Planar {2ni) ... Npi 

anar 

(2nfc) (19) 

This can be seen again in terms of contractions, because propagators joining different 
vertices must satisfy additional momentum constraints as opposed to the disjoint contrac- 
tions, therefore only the disjoint contributions survive in the large A limit. This amounts 
essentially to the cluster property. 
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Now recall that in the non-commutative case, the integral is given by the suitably 
normalized trace Hence the integrals above depend only on the eigenvalues of 0, and 
we obtain statements about the induced eigenvalue distribution. The above observables 
(|T!I|l can be written in the form 

{^Trr^)...{^Trr'=), (20) 

and completely determine the effective probability measure for the eigenvalues of 0. We 
can certainly write any such expectation value as 

{i^Trr'):.{^Trr')) = j d<p,...d<p,M<Pi, ^n) E E (21) 

i j 

for some measure which we would like to determine for large A. In order to 

absorb the infinities we introduce a scaling constant 

Uln(l + Aj), d = 2 



and write^ 
Then (jUj) gives 



= OiQip. (23) 

i^Tr^'} = I- (24) 

Now all expectation values (j^Tnp^^) are finite and have a well-defined limit A/" — > oo. We 
can then describe the eigenvalues of (p (in increasing order) by an eigenvalue distribution, 

ip{s) = ip,, s = j^, se[0,l]. (25) 

Then e.g. 

^ " '' ds f{^{s)) (26) 



H 

in the large M limit. The measure now becomes a measure /^[(/^(s)] on the space of 
(increasing) functions ^{s) : [0,1] i— > M. To find this measure /i[(/9(s)], we first note that 
the factorization property (fT^ implies that the measure is localized, i.e. 

(/(^,)) = /(^(.)) (27) 

for any function /, where Tp{s) is the (sharp and dominant) saddle-point or maximum of 
li[p]. This saddle-point ^(s) corresponds to a density of eigenvalues 

ds 

p(^) = — , / p{p)dp = l. (28) 
d^ J-oo 



^uq is not the wavefunction renormalization 
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The expectation value of the above observables is then given by 



(-^Tr/((p)) = J dppip) fip), (29) 

or (-^Try?") = Jq ds Tp^; for example, (/ dsip{s)'^) = J dsTp(sY = ^. We want to determine 
this saddle-point ^(s). This can be extracted from the above results: fTK|l implies that 

ds^isf'^) = j p{p)p'^dp = {^T Npianar{2n). (30) 

There is a unique eigenvalue distribution with these properties, given by the famous Wigner 
semi-circle law 

This means that the eigenvalues of are distributed correspondingly in the interval 

</)i e [-ao,ao]- (32) 
It is fun to verify (IHUj) explicitly for small n; indeed 

p{x)x^^dx = )l = - Npianar{2n). (33) 

_i i (2 + n) 4"- 

In general, this follows also from the basic properties of matrix models: Consider the 
Gaussian matrix model with action 



^0 = Urn) + ^Trcf 



(34) 



where = aoV? is a A/" x A/" matrix. Here fo{m) is some numerical function of m which will 
be determined below, and ao depends on m via c = c(m). 5*0 will find another interpretation 
in Section 1331 5*0 is known to reproduce precisely the eigenvalue distribution (j3T|l in the 
large M limit [19], and (jHOI) follows because again only planar diagrams contribute. Indeed, 
one finds again e.g. 

i^Tr^') = c^li^Tr^') = ^ «^ (35) 

etc.; we refer to the vast literature available on this subject, e.g. [19-21]. 

We conclude that if one is only interested in the eigenvalues, the free action 5*0 = 
/ d'^x^{di(f)di(j) + rn?(j?) can be replaced by the effective action Sq in Moreover, one 

can determine foiTn) such that the partition function is also recovered as 

Z = j ©06-^° = j V(Pe-^° = j r'0exp(-/o(m) - ^Tr^^); (36) 

this will be understood better in Section EISl To summarize, we noted that the observables 
(fTT?|l depend only on the eigenvalues; the factorization property implies that the dominating 
contribution comes from a well-defined eigenvalue distribution, which is via (jHn|) identified 
as the Wigner-distribution corresponding to a Gaussian matrix model. We can then write 
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down an effective Gaussian matrix model which reproduces all the expectation values for 
these observables. Note that the details of the propagator enter only through ao, and the 
basic result depends only on the degree of divergence. 

It is interesting to compare this with the commutative case, where the eigenvalues 
are replaced by the values of the field at each point, i.e. N'^ variables for a lattice 
regularization. Using again a cutoff A, we would have 

This could be reproduced by a Gaussian distribution = f d'^x(f>{x)'^. The crucial 

difference is that this would describe N'^ independent variables (which is false of course, but 
ok for these observables) with a very fiat distribution, while in the NC case we consider only 
TV" = N^^'^ eigenvalues which are collective and intrinsically non-local variables, governed 
by an effective action with a sharp potential^; note also the explicit factor A/" in (jH^)- This 
will become more explicit in the following sections, which do not apply to the commutative 
case. 

Let us try to interpret this result. It may be surprising that already the free scalar NCFT 
is apparently very different from the commutative case, even though they are supposed to be 
the same from the star-product point of view. The reason is that we are looking at statistical 
properties of the operator representation of the wavefunctions, which is related to their 
point-wise values only for low momenta where the functions are "almost-commutative". 
Therefore the non-classical properties of the high-energy modes are responsible for this 
property. This already indicates UV/IR mixing: the modes with k ^ are suppressed upon 
quantization because they tend to have the "wrong" eigenvalue distribution (in particular 
(f)k=o oc 1 has only one eigenvalue). 

The existence of a well-defined eigenvalue distribution and its matrix model description 
will generalize easily to the interacting case at least on a perturbative level, and is very 
plausible also in the non-perturbative domain. Furthermore it suggests new and practical 
insights to the coupling constants and their renormalization. This will be discussed next. 

3.2 Interactions 

If we include interactions of the form 

SU<P) = ^ / d'^xrix) = ^(27r^)<^/2yr0", 
n J n 

the results flH4|l . for the free case generalize at least on a perturbative level: consider 

J P0e-(^«+^'-) = J P0e'^o - ^{2neY^^Tr(l)'' + (37) 

Hence the expansion of e"*^™' in powers of Qn leads to additional terms of the form or 
()20|). which are again evaluated by the rule ()29p for the measure p of the free case^. We 

•^recall that the saddle-point approximation for the matrix models is good for the eigenvalues, but not 
for the matrix elements 

''This is valid as long as the coupling constants g„ are small enough so that the eigenvalue distribution 
remains unchanged. This should be enough to determine the perturbation series, where gn can be arbitrarily 
small. In particular, it suffices to determine mass renormalizations etc. as seen below. 



■'int 
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can write a similar formula to evaluate observables such as {jj'Tr(j)'^"')g in the presence of 
gn- Therefore the results of Section ITTl apply, and one can simply replace 5*0 by 5*0 in (jTf|l . 
We conclude that at least perturbatively, the eigenvalue sector of the NCFT with action 

S = j d''x^{d,(f)d,(j) + m^(f)^) + S,nt (38) 

for any polynomial interaction is described by the effective matrix model 

(39) 



^(0) = /o(m) + ^rr02 + S,rM) 



where ao = ao("^) is given by ()22|1 . in the large M limit. This is expected to be correct 
as long as the eigenvalue distribution is close enough to the free one (13111 . The partition 
function is similarly given by 

2J\f 



■'int 



V(f)exp{-fo{m) - ^T7-^Tr^' + S.nM))- (40) 



This defines by the usual matrix model technology an analytic function in the couplings 
Qn- In the later sections, this will allow us in particular to determine some renormalization 
properties of the potential is a very simple way. We will also explore non-perturbative im- 
plications such as phase-transitions, hoping that this will give at least qualitatively correct 
results. 



Scaling and relevant couplings. Let us try to estimate the impact of the interaction 
terms to the eigenvalue distribution. In the matrix-regularizations of NCFT used below, 
we will have M ~ A^'^/^ and A = 0{-\/N). Therefore the following scaling behavior holds 

0~ao~A^"\ N = K'^ (41) 

using (j^Hj) . We assume this scaling also in the interacting case in the weakly-coupled 
domain, and check for which Qn this is self-consistent. Then fl39|l essentially becomes 
S ~ Tri^Mip"^ + QnOiQip'^) =: Tr{MV{ip)), omitting constants of order 1 (including 6, which 
is not assumed to scale). Here if is of order one. The resulting eigenvalue distribution is 
governed by V{ip) [19], and it will remain near the Gaussian fixed point resp. Wigner's law 
provided the bare couplings satisfy Qn < A/'/ctg ~ A^, where 

S = d + n-'^ (42) 

is just engineering dimension of As usual, this means that relevant or marginal couplings 
with 5 > are "safe" and expected to be renormalizable, while irrelevant couplings with 
5 < must be fine-tuned and are expected to be non-renormalizable. 

In general, it will be safe to use as long as the shape of the resulting eigenvalue 
distribution Tp[s) is close to the Wigner law. For a qualitatively different shape, one should 
expect corrections to (jH^ . 

In particular, it is quite clear that turning on some small coupling (74 must be compen- 
sated by a suitable "mass renormalization" in order to preserve the shape of Tp{s). This 
allows to determine the mass renormalization, which will be explored in Section But 
before discussing these issues, let us look at the above results from a different perspective: 
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3.3 Angle-eigenvalue coordinates in field space 

One of the merits of the NC field theory is that it naturally suggests new coordinates in field 
space, which are very different from the usual "local" fields in the commutative case. 
This is particularly obvious using a regularization in terms of finite-dimensional matrices 
(p G Mat{M , C). Then there is a natural action of the SU {M) group — > U~^(f)U, which can 
be seen as NC version of the symplectomorphisms. Even though this is not a symmetry due 
to the kinetic term, it suggests to parametrize in terms of eigenvalues and "angles" , which 
are very non-local coordinates. This change of variables leads very naturally to the picture 
we found above. In particular, the non-trivial measure factor in the path integral due to 
the Jacobian makes the existence of a non-trivial eigenvalue distribution very plausible. 
We start from the simple fact that any hermitian matrix can be diagonalized, 

= [/-i(0,)[/ (43) 

where U G SU{Af) and (0j) is a diagonal matrix, with eigenvalues 0j. We can moreover 
assume that the eigenvalues of (0j) are ordered; then the matrix U is unique up to phase 
factors K = U{l)-'^~^, provided the eigenvalues 0j are non-degenerate. This leads to the 
following definition of the orbits 

(9(0) := {[/-i(0i)f/; U e SU{Af)} ^ SU{Af)/K (44) 

where K is the stabilizer group of {(pi). These are compact homogeneous spaces. Then the 
partition function can be written as 



P0 exp(-^(0))) = J d(t>,A\(t>,) J dUexp{-SiU-\(t>,)U)) 
= I rf0,rff/exp(5^1og 10, - 0,1 - S{U~\(Pi)U)) 

= j rf0,F(0,)exp(-(27r^)'^/2 ^ V^(0^) + ^ log |0, - 0^.|), (45) 

where dU is the integral over M ^ M unitary matrices; similar manipulations were also 
done in [9]. We introduced here the function 

F(0) := / dU eM-Sk.n{U~\cl>)U)) =: e'^^^), (46) 



which by definition depends only on the eigenvalues of 0. The last form is justified because 
F(0) is positive. It satisfies 

^(0 + c) = ^(0), ^(-0)=^(0). (47) 

Moreover, J^(0i) is analytic in the 0, because the space is compact, and invariant under 
exchange of the eigenvalues. We can therefore expect that it approaches some nice classical 
functional ^[0(s)] in the large M limit, where 0(s) is the function in one variable which is 
related to 0, as in (j2SI)- 



12 



We can now read off the induced action for the eigenvalues, 



i^j i 

The log - term in could also be absorbed by defining 

^(0,) = ^(0,)-^log|0,-0,|. (49) 

In particular, note that the log-term in (PH|) strongly suppresses degenerate eigenvalues, and 
this cannot be compensated by any analytic ^{4>i). Therefore the saddle-points of S'(0j) 
corresponds to some non-degenerate eigenvalue distribution. Furthermore the kinetic term 
(encoded in ^{(pi)) strongly suppresses jumps in this eigenvalue distribution^, therefore we 
expect it to approach some smooth function (j){s) in the limit M —>■ oo. Furthermore, we 
expect this eigenvalue distribution to have compact support after a suitable rescaling: 

0(s) = alp{s) (50) 

where a denotes the maximal eigenvalue. This typically happens for matrix models, and 
appears to be true also in this context as shown below. Indeed since J-'{4>) only depends 
on the eigenvalues, one can trivially interpret it as a function of any hermitean matrix 
= f/~^(0j)[/, and rewrite the partition function as an ordinary matrix model with formal 
U{M) symmetry, 

Z = j d0iexp(-^(0i)) = j P(/)exp(-J^(0) - (271^)^^/2 Tr\/(0)). (51) 

The last step is of course completely formal, and we are only allowed to determine ob- 
servables depending on the eigenvalues with this action. These are determined by the 
"effective action" S'(0i), since the degrees of freedom related to U are integrated out. This 
strongly suggests that much of the information about the quantum field theory, in par- 
ticular the phase transitions and the thermodynamic properties, are determined by S'(0j) 
and the resulting eigenvalue distribution in the large M limit. Note that this is essentially 
a one-dimensional problem, governed by the (unknown) functional jF[0(s)] and V . The 
advantage of this formulation is that it is very well suited to include interactions, and 
naturally extends to the non-perturbative domain. 

Now assume we know J^(0i); one can then look for the saddle-points of S[(j)i), determined 

by 

^ + {2^eY/' V'(0.) = (52) 

and ask if they are localized enough to dominate the observables. If the existence of a sharp 
eigenvalue distribution is established, the full path integral in ()45|) would be dominated by 
the integral over the corresponding SU{M) orbit 0{(f)i). This in turn should allow to recover 

^At strong coupling however, we will find a phase transition to a distributions with one gap, correspond- 
ing to some "striped" phase 
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also other properties such as correlation functions, by integrating over the corresponding 
which is compact. This will be discussed further in Section 
We can now relate this to the results of Section 13.11 All observables of the eigenvalues 
as considered there are determined by the "effective action" S{(f)i) above. The result of 
Sections Em and ini2l in particular the factorization property (jTHl), says that there is indeed 
a well-defined eigenvalue distribution in the non-commutative domain^, i.e. as long as 
^ 1. Comparing with the effective action for the free case ()34|) . we find 

2J\f 

H<P) = fo{m) + —Trcf" (53) 

which certainly reproduces all admissible observables for the potential V{(j)) = J |m^0^. 

This formula may appear strange, since ^{(f)) should of course be independent of V. 
The reason is that this relation (j53|) has been established only "on-shell", for eigenvalue 
distributions close to the Wigner law for the free case. It is not clear how well this works for 
large deviations from that case. Later we will use this form also for eigenvalue distributions 
which are quite different from the free one, where one should expect corrections to ()53|) . The 
appropriate way to use for some given eigenvalue distribution is therefore to determine 
m such that the corresponding free distribution matches best the one under consideration. 
It would be extremely interesting to know more about the functional J-'{4>)- 

The dominance of a given orbit 0{(j)i) in is clearly related to UV/IR mixing: 
naively, the action has a minimum at = const 1; however if the volume-factors from the 
path integral are taken into account (which happens at one loop), this zero- momentum 
state is actually highly suppressed, and the dominating field configurations have nontrivial 
position-dependence (i.e. momentum), due to the nontrivial eigenvalue distribution. Note 
that this argument is completely non-perturbative. Furthermore, depending on the form 
of the potential V{(f)) the dominating eigenvalue distribution may be connected or consist 
of disjoint pieces. These would clearly correspond to different phases. This picture will be 
made more quantitative below, and we will be able to identify the "striped" phase of [6]. 



3.3.1 Interpretation of the orbits (9(0i). 

Consider the reduced model pUj) for a given orbit O(0i) with fixed eigenvalue distribution, 
whose (free) energy is given by ^{<f))- Intuitively, we can interpret this model as follows: 
consider a classical fluid 0(x) on a compact space (due to the regularization e.g. on CP" 
or some torus) with prescribed "density" 

pip) = ^ j d''x6{<P{x)-p), (54) 

corresponding to the eigenvalue distribution. Then p{p) is essentially the density of eigen- 
values ()28j) . at least in the semi-classical limit. Note in particular that the action of the 
classical volume-preserving diffeomorphisms on 0(x) can be approximated by St/ (A/"), since 
any configuration with given p can be obtained using SU{N'). Entropy will favor mixing, 
but the "kinetic energy" suppresses mixing. In a small region of space, the global con- 
straint ()54|1 is quite irrelevant, and the fluid will behave like a fluid with the same action 

^Tliis is not true for conventional field theory with 6 — 0, even though the formulation of this section is 
still possible for e.g. the commutative limit of fuzzy spaces. 
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but without the constraint. However, if we fix the field on a large part of the volume, this 
must be compensated in the remaining space. This is clearly an IR effect, suppressing very 
large wavelengths. Therefore we expect that the theory behaves like an "ordinary" field 
theory in small enough regions of space and may hence describe ordinary local physics, 
however it is certainly different globally. This is quite interesting and encouraging for 
possible applications in elementary particle physics. 

In order to go beyond the computations in the following sections, one should therefore 
study the reduced models (pUjl in more detail, and see to what extent they approximate a 
scalar field theory. Note that quite generally if the scale a is increased, the configurations 
with short wavelength will be more strongly suppressed, leading a long correlation length; 
on the other hand for small a, the correlation length will be short. Therefore there should 
indeed exist some suitable scaling a{Af) which gives a macroscopic correlation length in 
the limit J\f — oo. The relation with the intuitive picture presented in Section and in 
particular (j5l|l can be seen best using coherent states or projectors (see e.g. [25,26,28]), in 
particular for the regularizations with CP^. It would be very interesting to combine these 
methods with the approach in the present paper. 

The "ground state" of on a given orbit 0{(j)i) is some very smooth function with 
the given density, which solves the e.o.m. 

[A</),0] = O. (55) 

This has many interesting solutions: One class is given by solutions of the free wave equa- 
tion A0 = c0. In particular, the (non-commutative) spherical harmonics with suitable 
eigenvalue density are solutions also on the above orbits 0(</)j). However there are other 
solutions, for example any solution of A(f) = f{(j)) for arbitrary / solves ()55|1 : in particular 
any diagonal matrix does (in the usual basis). A careful study of these issues should lead to 
a better understanding of (HH)), and hence to improvements of the simple results presented 
below. 

3.4 Relating the matrix model to physical observables 

Before analyzing further the matrix models (jH^ . we should try to relate them to the 
physically interesting quantities such as mass and coupling constants. Recall that the 
definition of mass on a non-commutative space is not obvious, since the Lorentz-invariance 
is generally broken. However on Euclidean Mg" with Sp{2n)- invariant 6ij one can define a 
correlation length in terms of correlation functions for 2 suitably localized Gaussian wave- 
packets^ (j){x) of size l/y/O, say, or e.g. coherent states for the fuzzy spaces CP^: 

(0(a;i) <P{X2)) = [ d<P,A\<P,) [ rff/e-^(^-'(*-)^)(0(a;i),[/-i(0,)[/)(0(x2),f/-^(0i)t/). 

(56) 

Here {(j){x),U ^{(j)i)U) oc Tr{(j){x)U ^{(f)i)U) denotes the inner product on the NC space. 
Note again that only the kinetic part of the action is non-trivial here, and we will assume 

^These "test-functions" can be moved transitively on the NC spaces using the translational symmetry 
which is unbroken. On Rg" with U(n)- invariant 9ij resp. their fuzzy versions CPJ^, the residual unbroken 
U(n) rotational symmetry is maximal and large enough to ensure that the correlation length is independent 
of the direction. 
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in the following that the full path integral is dominated by some orbit 0{(f)i). The mass 
can then be identified as the inverse correlation length, provided it is -C mg. Similarly, 
one could define the coupling constants in terms of correlation functions of e.g. 4 such 
wave-packets. In general, it is plausible that the interesting low-energy observables depend 
only on the eigenvalue distribution, and can be determined in principle by an integral over 
the orbit 0{(j)i). 

The question of renormalizability is then roughly whether one can scale the (finitely 
many) couplings Qn with the cutoff A ^ oo in such a way that the correlation length and 
all the other observables at low energy (or at the scale of non-commutativity set by 6) 
approach a well-defined limit. To answer this question of course requires control over all 
these correlation functions. In this paper, we try to proceed as much as possible without 
resorting to perturbation theory. In view of the above results on the eigenvalue distribution, 
it seems very plausible that the scaling of the bare couplings m and Qn must be determined 
such that the shape of the dominating eigenvalue distribution, i. e. the normalized function 
^(s) fixed as A ^ oo. These scalings should be easily accessible with matrix model 

techniques if we know the matrix model, in particular J-'{(f)) . Moreover, it seems plausible 
that this should even be sufficient to guarantee "renormalizability" in the above sense, 
provided the field theory can indeed be reduced to the orbit O(0i) as discussed in Section 



To make explicit computations, we have to use pTIjl resp. for the time being, i.e. 
we have to require that ^(s) respectively its related eigenvalue density (j2Hl) is close enough 
to the Wigner law ()31|). Then the "physical" mass mpt resp. correlation length can be 
identified without computing such correlation functions: it should be given by the bare 
mass of the corresponding free theory with the same maximal eigenvalue ao{mR), for the 
same cutoff A. That is, mn is determined by 



where a ()5()j) will depend on the couplings of the interacting matrix model. This will be 
elaborated in Section |3] We will apply this prescription (jKTjl even if the function ^(s) is 
not close to the free one in this paper, hoping that it is mainly the "size" of the eigen- 
value distribution which determines the correlation length. This is plausible in view of the 
classical picture discussed in Section 13.3.11 

We want to point out again the following consequence of this picture: (j56|l is strongly 
suppressed for zero momentum 0o 1, since {(f)o,U~^{(j)i)U) = for the dominating 
eigenvalue-distribution. On the other hand for non-zero momentum, the eigenvalues are 
non-degenerate, and the above inner product is non-zero. This suppression of zero mo- 
mentum is clearly related to UV/IR mixing, and suggests that indeed only localized wave- 
packets should be used. 

4 Example: the 0^ model 

Consider now the model 



a = ao{mji) 



(57) 



S 
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d^x 09,09,0 + ^m^02 + Q Am202 + |04^ ^ (58) 

where we have introduced a "renormahzed" mass m\ = — Am^. Its eigenvalue sector 
according to the above results is described by 

~S = foimn) + ^^^^^ Tr02 + {2n9Y/^Tr (^^Am'^' + |0^^ . (59) 

While m/j is arbitrary in principle, we will adjust it in order to minimize the expected 

errors due to our only partial knowledge of J^(0i). Note that in the regularization using 
fuzzy CP", we will find (ITHi|l 

A=V^ (60) 



where N is related to Af via ()145|) . and a similar result A = using NC tori pi8|) 

where J\f = N'^^'^ . Therefore 

^ = 0(iVe('^-2)/2) = OiA'e''/^) (61) 
"o 

up to log- corrections. 

To solve this model, we first rescale as 

= <yg(p (62) 

such that the saddle-point solution for (p will have an eigenvalue distribution with range^ 
[—1,1]. Rewriting the matrix model ()59|1 as 

S = foimn) + MTr {^^^ + ^-^TvA . (63) 



we have 

(2i,0f^gal = A/-J', (65) 
= m\ + Am^. (66) 

The model ()63p is well-known and can be solved with standard methods from random matrix 
theory, see e.g. [19,21]. It is again governed by a single saddle-point in the eigenvalue sector, 
with eigenvalue distribution ^^(s) resp. Pg{p)- Note that now m'^ < is allowed. Assuming 
that g' > 0, there are 2 cases corresponding to distinct phases of the model which we are 
discussed below. 

It is interesting to note that even small g' < is admissible for this matrix model as 
long as m'^ > 0. This indicates analyticity in g', so that perturbation theory does make 
sense in the weakly coupled phase. 

^in the disordered phase discussed at present, see Section H!T1 
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4.1 Phase 1 ("single-cut"): m'^ > 0, or m'^ < with f-7 < 1 



In this case the eigenvalue density is given by [19] 



pM = ^iaV + ^ + m")VT^'. (67) 

We have imposed that the range of be [—1, 1] as explained above, which using standard 
formulas [19] implies 

1 = 7^(-m'' + ^m'^ + 12g') (68) 
^9 

i.e. 

g' = ^(A-m'') (69) 

Note that the matrix model ()63|) has 2 independent parameters, which can be chosen either 
as m'^ and g', or e.g. ag and g'. We have chosen to work with ag, therefore m'^ and g' are 
not independent. In particular, note that ^ < 1 for g' < 16 due to (jEHI), therefore this 
phase will be the "weakly-coupled" phase in Sectional 
Inserting m'^ = 4 — | (7' in ()64j) and using ()65|) gives 



a. 



+ {2ne f'al^Am' = 2 - - {27:6^/' J g. (70) 



_g_ 




So far, m/j was arbitrary, and it shouldn't matter for fixed "bare parameters" m and g. 
However, the relation (jHU)) is expected to be good only if the eigenvalue distribution is close 
to the free one corresponding to m^j. We should therefore choose mpt accordingly, and a 
good choice is 

aoirriR) = Ug (71) 

as in ()57|) . This guarantees that the eigenvalue distribution of the interacting model (jU^ 
has the same range as the free one obtained from mn, hence it is "close". Then (f7n|) 
simplifies further, and using Am^ = — m|. we have 

{2n9f'ga',{mn) = Afg' (72) 
+ ^9 alimR) = m\. (73) 

These are 2 equations in 4 unknowns (the function ao{mR) being known ()22j) ). and we can 
basically choose any 2 of them to parametrize the model. 

4.2 Phase 2 ("2 cuts"): m'^ < with ^ > 1 

At 4(7' = m'^ with m'^ < the eigenvalue density breaks up into 2 disjoint peaks, which 
are concentrated around the 2 minima. For Ag' < m'^ these peaks have finite distance, 
and moreover spontaneous symmetry breaking of Trcj) occurs with non-zero "occupation" 
in both peaks [22] . This clearly describes a different phase of the model. The semi-classical 
interpretation in the picture of Section l!^.3. II would be a 2-fluid model, with a large surface 
energy at the contact surface due to the kinetic term. It is very plausible that this will 
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develop some kind of (generalized) striped pattern upon mixing, and we conjecture this 
to be the "striped" phase of [6], or equivalently the "matrix phase" of [9]. Note that the 
assumption in [9] of 2 delta-like peaks is qualitatively very close, but strictly ruled out by 
the log-terms in PH|) . 

It is even possible to consider negative g' in this model [19], reflecting the analyticity in 
g' . If g' becomes too negative there is another phase transition, which will not be considered 
any further here. 

5 Weak coupling and renormalization 

We first consider the weak coupling regime, and try to get some insights into the renor- 
malization of the 0^ model. In particular, one would like to understand better the relation 
of the above approach with the usual concepts of running coupling constants etc. The 
natural scaling parameter here is the size of the matrices M ~ A^'^/^, which is related to 
the UV-cutoff A through (jHIll) for the regularization using fuzzy CP". We can therefore 
ask how the bare parameters m?,g etc. must be scaled with A such that the "low-energy 
physics" remains invariant, and we could define the corresponding beta-functions. 

In the weak coupling regime, the shape of the eigenvalue distribution Tp resp. Pg 
(jU7j) should certainly be kept close to the free one, given by the Wigner law. As discussed 
in Section 13^ the correlation length resp. "renormalized mass" itlr can then be obtained 
simply by matching the size of the eigenvalue spectrum ag with the one for the free case 
with the same cutoff, 

ag = aoirriR) (74) 

see ()57|1 . Then mn measures the correlation length of the free model which best fits the 
eigenvalue distribution. This should be very reasonable, since we are looking at a weak 
coupling expansion where g' can be taken arbitrarily small. 

The corresponding running of the original parameters m and g depends on the dimen- 
sion, and will be worked out below. For the mass, this is very easy in our approach. 
Unfortunately, it is not obvious how to relate the coupling g resp. g' to e.g. 4-point func- 
tions or scattering processes; this could be supplemented perturbatively by a conventional 
RG computation, cp. [11,12]. Also, note that there is no wavefunction-renormalization in 
this approach, as there are only 2 coupling constants in the matrix model (jHSJ- On the 
other hand the computations below are much simpler than any standard field-theoretic 
computations, and provide non-trivial information on the model without having to resort 
to perturbation theory. 

5.1 2 dimensions 

To find the scaling of and g we can use (f7^ and (fTHj) . where ag = ao(mji) is already 
built in. The running of the mass can be read off from ()7H|1 . 

3 3 A^ 

= ml- -galirriB) = ^R-^9 ln(l + ^) (75) 

Recall that rriR is the "physical", renormalized mass which is supposed to be finite while 
A ^ cxD. We see explicitly the expected logarithmic divergence of the bare mass. Note that 
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this is an all-order result in g, despite appearance. 

It is interesting to compare this with a conventional one-loop calculation. Noting that 
our interaction term is J (which is more natural in the matrix model context) rather 
than J ^0^, one would find 

c 2 g r (fk 1 3 , , A2 

5m,^^, = 12- j = j^aHl + ^) (76) 

from the 8 planar plus 4 non-planar contractions. It is interesting to see that this one-loop 
computation agrees precisely with our procedure of matching the eigenvalue distribution 
with the free one, even though we used the non-perturbative matrix model results (jU7j) 
ff. One might have expected that only the planar diagrams contribute to the mass renor- 
malization; however since the "mass" is obtained as the p ^ limit of the IPI 2-point 
functions (for fixed A), the planar and non-planar diagrams coincide and both contribute. 

Note that if we use another regularization such as NC tori, we should use the modified 
propagator ()128|) rather than the sharp cutoff in This would lead to a somewhat 

modified formula for q;o(?77.r), however the corresponding mass renormalization ()75|) would 
still coincide with the one-loop result since the same ao{mji) enters in both calculations. 

Next consider the bare coupling g. Since the matrix model provides no simple relation 
between g and the 4-point function, the latter would have to be obtained e.g. by a per- 
turbative calculation as usual. However since g is not expected to run in 2 dimensions, we 
simply interpret g as "physical" coupling constant, fixing the scale by rrig. The relation 
with the matrix model coupling g' is given by (f7^ . which for large A and using (jUUI) is 

In particular if we keep g finite, we see that g' ^ very rapidly as A ^ oo. This means 
that the eigenvalue density pg (jHTj) is very close to the free one given by Wigner's law, 
and we expect that the above relations are reliable in this weak-coupling phase of 0^ in 2 
dimensions^ . 

5.2 4 dimensions 

Using the same procedure, (f7H|l gives now the expected quadratic divergence of the bare 
mass 

lovr^ \ A^ mj^ J 

This is again an all-order result, and agrees again perfectly with a conventional one- loop 
computation (see e.g. (3.11) of [5] after replacing g — >• gj^. However the bare coupling g is 
now also expected to run, being logarithmically divergent at one loop. Recall also that we 
cannot see any wavefunction-renormalization as there are only 2 parameters in the matrix 
model. The relation between g and the matrix-parameter g' is now given by 

= - 1-2-fln— g^-g 79 



^Recall however that the non-planar diagrams are suppressed only logarithmically (|14() in 2D, which is 
much weaker than in 4D 
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for A ^ mji, where fl60p and A/" = ()145p has been used. This means that (keeping 

in mind some possibly logarithmic renormalization of g) our procedure to derive ()78|) is 
justified in a perturbative sense, since we can make g' arbitrarily small so that again pg is 
very close to Wigner's law. 

It is quite remarkable that according to the above results, the mass is renormalized 
only at first order in the bare coupling g, even though the derivation is non-perturbative. 
While this is expected for d = 2, this seems very surprising for d = 4; it would indicate 
that NC (j)'^ in 4 dimensions is much better behaved than in the commutative case. This 
can be traced back to our procedure of determining the mass m/j through a = ao{mfj). 
While an exact calculation of the correlation length on a given orbit 0{(f)) would certainly 
modify this result, the basic picture seems clear and simple. Unfortunately this analysis 
does not give us the relation between g and the "physical" coupling as obtained by the 
4-point function, which probably requires further renormalization. 

It is tempting here to conjecture that fixing g' and m' or equivalently the eigenvalue dis- 
tribution ^ in a suitable way suffices to define a non-trivial NC field theory in 4 dimensions. 
We will find further evidence for this in Section Ifi.ll 

We conclude that while this approach certainly needs further work and thus far provides 
only a partial window into NCFT, it suggests a very simple and compelling approach to 
renormalization in the NC case. This may help to overcome the difficulties at higher order 
found in [5]. 

5.3 Higher dimensions 

It is also illuminating to try to generalize the above considerations to dimensions higher 
than 4, where the commutative 0^ model is no longer renormalizable. The relations (jUBj) 
and (fT2|l are still valid, and give 

m^ = ml-g 0{A''-^) (80) 

and 

g' = 0{A'-') g. (81) 

Clearly our matrix model (jfjHj) makes sense only for g' = 0(1), which would require g = 
0(A^^^). This is just the canonical dimension of g. If we again assume that g' and hence 
the shape of the eigenvalue distribution should be independent of A, we see that g must be 
fine-tuned to at least 0(A^~°'). This is in accord with the expected non-renormalizability. 

6 The phase transition. 

Now consider the phase transition of the matrix model (jHSl); which is known [21,23] to have 
a phase transition between the 1-cut and the 2-cut phase at 

m'^ = Ag'. (82) 

This is expected to be the phase transition between the disordered and striped phase of 
Gubser and Sondhi. Combining with this gives 

g'=^(A-m'') = ^m'' (83) 
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which has 2 solutions for m'^; only the negative one 

= -8 (84) 

is relevant here and marks the phase transition. The corresponding coupling is^'' 

g' = 16. (85) 

Note that there is still a free parameter in the model, which we can take to be rriR resp. a. 
Plugging this in (|7^ gives 

{2neY'''gal{mR) = 16AA, (86) 
^'^ + ^9 (^{mR) = m%, (87) 
and solving the first equation for we obtain 



'^" + ^^J0fn=<^9l (88) 

where m|j is a function of g via ()86j) . This defines the critical line. The corresponding 
phase transition is third order [23] in the variables g' and m' . This is in contrast with [6-8], 
who argued for a first-order phase transition. 

Note that m'^ < means that we are quite far from the perturbative domain, and the 
eigenvalue distribution is significantly changed from Wigner's law pip. In this regime the 
replacement ()34|1 for the kinetic term has not been tested, and we cannot expect the results 
to be exact. Furthermore, we expect the "bare" parameters m,g to have some non-trivial 
scaling in A/" ~ N'^^'^ on the critical line corresponding to the RG flow, which is not obvious. 
In any case, it is quite reasonable that this description is at least qualitatively correct. To 
proceed, we have to discuss the dimensions separately, starting with the most interesting 
case of 4 dimensions. 



6.1 4 dimensions 

We consider first the regularization using fuzzy CP^ discussed in Appendix A, which cor- 
responds to a sharp cutoff. In this case we have M = N'^/2, and = O(A^) in d = 4 (j221). 
By looking at the equations (jHEI), (|H7jl and 

A= ^ f (89) 

which holds for this regularization of M.g, it is quite obvious that there should be solutions 
which scale as 

oc A^ oc A^, g fixed. (90) 

Hence we expect a phase transition at finite coupling g and the standard quadratic running 
of the mass. To find a closed equation for the critical line, we have to use (j^^ for ao = 

^J^^ = l^!Ain(l + ^) (91) 
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Recall that < 16 is the weakly coupled single-cut phase as pointed out in section ITTI 
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This cannot be solved explicitly for m^, but we can use ()88p for which together with 
M = and ^ gives 

w? ?) [q 171% , , 

Since m|j > 0, this makes sense only for 

vn^ 3 [g 

Plugging this in and using 



.o > -— a/-- (93) 
A2 47r V 2 ^ ' 



9/ X ^ /2 

=^Vr^ '''' 



we get 



This is indeed consistent with the scaling (jHSI). Since the rhs is G (0, 1], this has a solution 



only^^ for 27r^/| < 1, i.e 



11 f_ /2 ^ 

2 



^ > ^7c = Stt^. (96) 

In terms of the dimensionless parameter = m^/A^, the corresponding critical mass is 
ml = — §• Expanding g = gc + 6g, w? = ml + dw? the critical line ()95j) is given by 

5(; = _^5m2 + ... (97) 

for small variations. This is plotted in Figure ^ 

Let us try to assess the validity of this result. The basic arguments underlying this result 
are expected to be good in 4 dimensions as long as m|j <^ A^ (unlike in 2 dimensions, see 
below and Appendix B). This is satisfied here since as g ^ gc from above, (j^HI) implies that 

using (jn2I), and in particular mR = at g = gc- Therefore the replacement of the free action 
with the matrix model (jH^ is essentially justified, apart from the modified eigenvalue 
distribution (p{s). It would be very interesting to estimate the effects of this modified 
^(s) more rigorously, e.g. by estimating the integrals over the compact orbits 0{(f)i). 
The existence of gc can also be understood simply by noting that the critical line is 
characterized by a specific eigenvalue distribution which is different from Wigner's law, 
which holds for g = 0. Therefore the critical line^^ cannot end at g = 0. While the 
existence of a critical point should be a sound prediction, the relation between and g 

^^the rhs of (|95|l is simply c(A, mji) as defined in 0, which will very generally satisfy such bounds. 

""^^note that rn^ cx A'^ here, as in the weakly-coupled regime. This is different in 2 dimensions, where the 
critical line is in a different scaling regime from the weakly-coupled case, and no such conclusion can be 
drawn. 
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Figure 1: The critical line in 4D with critical point 



on the critical hne cannot be expected to be precise. If we would use the propagator for 
the lattice regularization ()128|1 . the details of (jHSI) would change but the results would be 
quahtatively the same. In particular the critical coupling could be evaluated using ()132j] . 

The existence of a critical point terminating the critical line at Qc ^ is certainly 
intriguing. Since the critical line will be stable under the RG flow, its endpoint should 
correspond to a non-trivial RG fixed point for the non-commutative 0^ theory, and the 
correlation length is expected to diverge^^. That limiting model should be non-trivial, 
since the eigenvalue distribution Tp[s) is different from the free case. It is furthermore 
plausible following the discussion in [6] that it admits a "continuum" interpretation as a 
renormalizable NC field theory. The physical content of such a fixed point would require 
to study e.g. 4-point functions and the running of the coupling, which is beyond the scope 
of this work. Note however that if we assume that g becomes larger with increasing cutoff 
as in the commutative case, this would mean that the low-energy coupling corresponding 
to Qc is small, and hence in a physically interesting regime. 

These results and in particular the existence of a phase transition are roughly consistent 
with the results of [6-8], but not precisely; for example [6-8] argue for a first-order phase 
transition using self-consistent Hartree-Fock approximation resp. a one-loop RG analysis, 
while we find a higher-order transition using a non-perturbative matrix model result. The 
approaches are thus quite orthogonal: in the present approach the interaction is treated 
exactly while the kinetic term is approximated, whereas [6-8] do the opposite. The fact 
that one always finds a phase-transition is quite encouraging. It is also remarkable that 6 
does not enter our result recall however that we always assume ^ 1 and ()154j) . In 
particular, one would expect a standard Ising transition to a uniform symmetry breaking 
state at small 6A'^; however our approach is not valid in that regime, and we should therefore 
not expect to see this transition. 

In dimensions higher than 4, ()86p together with = 0(A'^^^) = 0{N'^^'^^^) implies that 

^•^we find indeed mji = at gc, which should however not necessarily be taken at face value 
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g = 0{N 2+2) must be fine-tuned in order to stay in the weakly-coupled phase. This is 
consistent with the expected non-renormalizability, and we will not pursue this any further. 

6.2 2 dimensions 

In 2 dimensions our approach is more delicate as the planar diagrams are only logarithmi- 
cally divergent. To be safe we should allow only finite m^j while A —>■ oo, see Appendix B. 
Using and (IHUj) for = ao{mji) and Af = N, we have 

ml = 5 = = . 99 

^ e™o-l 2>/Ma ^ ^ 

e V 9 — 1 

Plugging this in ()88|1 and using the relation = 2^ for the regularization using the fuzzy 
sphere gives 

2,3/^, A^ 



m^ + - J- A = = = mi. (100) 
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This would be consistent with a scaling 

m^ ~ A^ A^ (101) 

on this critical line^^. However, the basic assumptions of Sections 13. II and Appendix B are 
no longer valid in this case, since then m\ oc A^. Our approach should be most reliable 
if mR remains finite, which implies through that g ~ 2(j^)^ and therefore pOOj) 
~ ~f ifx' Si^^^ cannot solve (jlOOj) in closed form, let us assume that g is smaller 
but not much smaller that g ~ 2(j^)^, so that we can neglect the term on the rhs of ()1()()|1 . 
Then the critical line should be given approximately by 

2 3 [g 



™^ = --^|A. (102) 

or a slightly modified formula using the fuzzy torus regularization (due to the different 
propagator). Unfortunately we cannot compare this with the numerical results of [10], who 
consider a different scaling. However we can compare it to some extent with the numerical 
results of [9] on the fuzzy sphere, and find reasonable agreement. This will be done in the 
following section. 

If g scales like fllOip . one can use a simpler argument due to [9] neglecting the kinetic 
term altogether, and replace the action by the pure potential model TrV{(f))\ hence this 
phase was denoted as "matrix phase" in [9]. Note that indeed the scaling (jlOlj) is appro- 
priate for the matrix model (jKH|) without the kinetic term. This amounts to identifying 



m'2 = 27r^— , g' = 2n9^, (103) 



which would predict a phase transition at {2n^^y = Svr^, i.e 



AT 7 ^" AT! 

m^ = --j2^A. (104) 

TT 



^''note that both and g have dimension of mass'^; nevertheless, this is a rather strange scahng in 2 
dimensions. This phase transition is very far from the weak couphng phase 
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6.2.1 The fuzzy sphere 

For the case of the fuzzy sphere, we consider the action 



S = l^Tr(0A0 + + A0^) (105) 



using the (redundant) parameters r, X, R following [9]. The eigenvalues of A = ^ are 
/(/ + 1)/-R^ = p^, with cutoff A = Pmax = N/R. Comparing with (|^. the above parameters 
are related to the coupling constants m and g in ()58|) via 

^ = ^, 2r = m\ 4:X = g. (106) 



It turns out that for the fuzzy sphere p4H) agrees precisely with the result for 



d2 



a^(m^) = iln(l + ^). (107) 



We can therefore apply ()102|) . assuming that A <^ (i^)^! particular we can fix A = 1 as 
in [9]. Then the critical line is 

r 3 1 1 

— = ^ - ^ -0.846 -. (108) 



Note that i? is a free parameter corresponding to 6, and we should take at least R = 0(v A^). 

Let us compare his with Martin's disordered-matrix transition [9] : he finds numerically 
a phase transition for (see eq. (46) in [9]) 

^ ^ -0.56 ^ (109) 

for large R. This is reasonably close to our analytical result. Recall that we cannot expect 
our prediction ()108|1 to be exact since the eigenvalue distribution is quite far from Wigner's 
law, and moreover the arguments in Section 2 are weaker in 2 dimensions compared to 
4 dimensions since the crucial divergences are only logarithmic. Furthermore g strongly 
violates g ~ 2(j^)^ and therefore mji —> 0, so that the replacement of the kinetic term 
with (jH^ is not justified here. The pure matrix model discussed above p03|l would give 

r V2 1 1 , , 

Hence the numerical results of [9] (taken for A^ ^ 30) appear to be in between the pure 
potential model and our approach, and our treatment apparently overestimates the kinetic 
term. This is not too surprising in 2 dimensions in view of the above remarks. 



7 Discussion and outlook 

We presented a simple non-perturbative approach to scalar field theory on Euclidean non- 
commutative spaces, based on certain matrix regularizations of Mg". Starting from a repre- 
sentation of the field (p = U{(f)i)U~^ in terms of eigenvalues (pi and "angles" f/, we observe 
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that the different l:)ehavior of planar and non-planar diagrams due to UV/IR mixing implies 
a particular eigenvalue density distribution, which can be reproduced by a simple matrix 
model. This is shown starting with the case of free fields, which can be described by a 
Gaussian matrix model. Interactions of the form TrV{(l)) can then be included very easily, 
and modify the eigenvalue distribution. This leads to a picture where the basic properties 
of the QFT such as correlation length and renormalization are related to its eigenvalue 
density p{s), through a "constrained field theory" with compact configuration space 0{4>). 
It also makes the existence of well-behaved scahng limits of NCFT i.e. renormalizabihty 
very plausible. We conclude that the eigenvalue sector of non-commutative scalar field 
theories is goverened by an ordinary matrix model, which provides a simple and intuitive 
window into the non-perturbative domain. 

For weak coupling, this approach provides a new way of computing the renormalization 
of the potential; in particular, using a very simple approximation we found an expression 
for the mass renormalization, which coincides with the conventional one-loop calculation. 
Furthermore, we found a phase transition at strong coupling in the 0^ model in both 2 
and 4 dimensions, which is identified with the striped or matrix phase of [6,9]. This is 
particularly interesting in the 4-dimensional case since the critical line then terminate at a 
non-trivial point with ^ 0, which is interpreted as an interacting RG fixed-point. The 
existence of this fixed point can be understood simply by noting that the critical line is 
characterized by a specific eigenvahic distribution which is different from Wigner's law, 
which holds for g = 0. Therefore the critical line cannot end at g = 0. All this suggests 
that such NC field theories may in fact be more acccesible to analytical tools than their 
commutative counterparts. 

Perhaps the main shortcoming of this approach is the lack of a precise relation between 
the eigenvalue distribution and the relevant physical parameters, such as correlation length 
resp. mass and coupling strength. It is quite clear that the leading parameter is the size a of 
the maximal eigenvalue, which has been used in this paper to extract physical information. 
To go beyond this approximation may be difficult, and may require e.g. perturbative 
methods. 

There are many other questions and gaps which should be addressed in future work. 
For example, the assumption that 6ij is non-degenerate (or special) is quite clearly not 
essential, and a similar approach should also work in odd dimension. Furthermore, a more 
elaborate analysis of the renormalization in the weakly-coupled regime should be attempted. 
Another interesting question concerns the relation of this approach with the results of [11] 
on a modified 0^ model; to address this issue, the above analysis should be repeated with 
a suitably modified propagator. This will be done elsewhere. 

Perhaps the most interesting perspective is the possibility that careful estimates of the 
contribution of the kinetic term on the orbits 0{(j)) (for eigenvalue distributions different 
from Wigner's law) should allow to rigorously justify the above picture also in the non- 
perturbative domain, and in particular the existence of the critical point. One may in 
particular try to establish renormalizabihty in this way. This should be facilitated by the 
fact that 0{(f)) are compact spaces. 

Finally, it would of course be extremely interesting to compare all this with numerical 
results in 4 dimensions, which are not available at this time. 
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Appendix A: Regularizations of 

The fuzzy torus and x Tjj 



)2n 



A particularly simple regularization of which works in any (even) dimensions was pro- 
posed in [18], which we review for convenience. Consider a (toroidal) lattice with lattice 
constant a and N sites in each dimension. We denote its size with 



L = Na. 



;iii) 



Since we are on a torus, one should not use the unbounded operators Xj. Instead consider 
the unitary generators 

(112) 



zN 



1. 



The commutation relations = iQij then become 

47r2 



ZiZj = exp 



L2 



i^ij 1 ZjZi. 



Rather than going through the most general case, we simply consider 

^ij ^ Qij 

and work out the 4-dimensional case, where 



(113) 



(114) 



Q 



1-j 





1 





o\ 


-1 




















1 


\0 





-1 


V 



;ii5) 



The generalization to any even dimension is obvious. Periodicity implies a quantization of 
a resp. 9 as 

^-e, (116) 



assuming that N is odd. The physical momentum is 

27rnj , TV TV . 
h = —r^ G— ,- = 
aJM a a 



-A, A), 



with UV-cutoff at 



A 



TT 



a V 9 

This is in qualitative agreement with the scaling ()154|1 obtained using fuzzy CP". 



;ii7) 



:ii8) 
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In order to write down the action for a scalar field, we also need partial derivatives or 
shift operators, 

Dj := e"^^ DjZiD] = e^^^'^^^/^Zi . (119) 
In our case, they can be realized as 



(120) 



etc. A solution of ()113|1 and Zj^ = 1 in 2 dimensions is given by the unitary "clock and 
shift" operators (recall that is odd) 



/O 1 
1 



\ 



1 
/ 



n 



\ 



^2(i-Ki/N) 



'1211 



which extends to arbitrary even dimensions by taking tensor products. Hence the field 
is a hermitean N'^/'^ x N'^^'^ matrix. The integral is again defined as 



/ = {2'ndY''Tr{f). 



(122) 



One can then define the "plane waves" 



d 



/Af^ 



i=l 



which satisfy 



Tr{(t>i(t>ri') = 6fin', (pi = <P-n 



(123) 



(124) 



for Hi G [— (A^ — l)/2, (A^ + l)/2]. They form a basis of the "space of functions 
Mat{N'^'^,C). Using 

Di(j)fiD\ = exp{2iTini/N) cj)^ 



one can write down the discretized lattice- version of ((H), 



]0) + ^0^ + f0^ 



(125) 



(126) 



For hermitean = '^Pn(pn with = pt-, the kinetic term becomes 

i 



(Jj 



(127) 
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The propagator is therefore 

= hi. ^2,(i-ooKM) + ".^' 

and the phase factor for the nonplanar diagrams is obtained from 

(Pk^Pk' = exp(-i^ Yl ^iQiiK) ^k'<Pk- (129) 

i<j 

ckq on the fuzzy torus 

Due to the different behavior of the propagator for large momenta, Oq on the fuzzy tori will 
be somewhat different from regularizations using a sharp cutoff, such as on fuzzy CP". In 
2 dimensions, we should compute more carefully 



<fx<t>\x)) = V / 

J —Tr/a 



V 



la (2vr)2 2E,(l-cos(p,a)) + mV 
" 1 



(27r)'^ - cos(ri)) + m2aV2 



(130) 



and in 4 dimensions 



-.la (2^ 2E,(l-cos(p,a)) + m2a2 



y .1 r d^r 1 



Jo (27r)4 ^.(l-cos(rO) + m2aV2 



(131) 



for the above regularization. In particular, for m = one has 

2/ ^ 1, . 4: . r d'^r 1 0.31 ,2 / X 

m < alim = = — / -—7 = -— A^ 132 

numerically. This is needed e.g. to determine the critical coupling Qc for the 0^ model in 4 
dimensions. 

The fuzzy sphere 

The algebra S^f of functions on the fuzzy sphere [16] is the finite algebra generated by 
Hermitian operators Xi — (xi, X2, xz) satisfying the defining relations 

[xj, Xj] = iANeijkXk, (133) 
xj + xl + xl^ (134) 

where it! is an arbitrary radius. The noncommutativity parameter Ajv is of dimension 
length, and is quantized by 



R -1 
30 



This can be easily understood: ()133|) is simply the Lie algebra su{2), whose irreducible 
representation have dimension A^. The Casimir of the A^-dimensional representation is 
quantized, and related to by ()134p and p35|) . Thus the fuzzy sphere is characterized by 
its radius R and the "noncommutativity parameters" or Aat. The algebra of "functions" 
Sff is simply the algebra Mat{N) of NxN matrices. It is covariant under the adjoint action 
of SU{2), under which it decomposes into the irreducible representations with dimensions 
(1) © (3) © (5) © ... © {2N — 1). The integral of a function / G 5*^ over the fuzzy sphere is 

0(a;) = ^Tr[0(x)], (136) 

which agrees with the integral on 5*^ in the large limit and is invariant under the SU (2) 
rotations. The dimensionless coordinates Aj = Xi/Aj^ generate the rotation operators Jf. 

J^f=[K,f]. (137) 

One can now easily write down actions for scalar fields, such as 

S = ^Tr + ^m'<p' + = 2n9 Tr Q^A^ + Im^' + \g<P'^ (138) 

were is a Hermitian matrix, and Acj) = -^JiJiCp is the Laplace operator. The last form is 
obtained by defining 6 through 

^ NO , , 

R' = -Y- (139) 

There are 2 obvious N ^ oo limits: 1) the conventional, commutative 5*^ limit keeping R 
fixed, and 2) the limit of the NC plane Rg, which can be obtained by keeping 6 constant: 
then the tangential coordinates Xi^2 satisfy at the north pole the commutation relations 



[xi,X2] = 'i'-^xs = i^^jR^ -xl-xl^ iO. (140) 

Note that 9^^ now determines the basic (NC) scale of the NC field theory, replacing the 
radius R. 

If one considers non-planar loops, the oscillatory behavior due to the 6j symbols sets 
in for angular momenta l"^ N and was studied in [24]. Note that this corresponds to 
= 0{6~^) in the NC case, but to p ^ oo in the commutative case. Therefore the non- 
oscillatory domain is divergent in the commutative limit and prevents a matrix behavior. 
This "low-energy" sector is suppressed in the NC case if 6 is kept finite. 

for the fuzzy sphere 

The arguments of Section 13.11 for the eigenvalue distribution go through here as well pro- 
vided the oscillating behavior of the non-planar diagrams is strong enough. Even though 
there are nontrivial effects even for finite R resp. 9 (xl/N [24], they are sufficiently strong 
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for our purpose only if 9 is kept finite. To find the appropriate aQ{m), consider for = 



^/(/ + l)/i?2 + m2 7q x2 + m2 

A^2 V A2 
R^ ln(l + — — ) = — ln(l + — ) (141) 



rn?R'^ An m 

in the large N limit, using V = AnR^ and x = ^-^^ and A = ^. This agrees precisely with 
(|TT|l . so that the same ao{m) as in (j221) can be used. 



Fuzzy CP" 

The construction of fuzzy CP" [27-30] is analogous to the case of fuzzy 5*2 = CP^. Since 
CP"" is an adjoint orbit of SU{n + l), it is a compact symplectic space and can be quantized 
in terms of finite matrix algebras Hom{VN) where Vat are suitable representations of su{n + 
1). To identify the correct representations Vn of su{n + 1), we must match the space of 
harmonics on classical CP" 

oo 

C°°(CP") = © V^p,o,...,o,p). (142) 
with the decomposition of Hom{VN)- It is easy to show that indeed 

N 

Hom{VM) = Vm®V;^ ® ^(p,o,....,o,p) (143) 

for 

Vn ■= V(Ar^o,...,0)- 

Here V(/j,...,/„) denotes the highest weight irrep of su{n+l) with highest weight /iAi+...+/„A„ 
where A^ are the fundamental weights. One can therefore define the algebra of functions 
on the fuzzy projective space by 

CP^ := HomciVN) = Mat{M, C) (144) 

with 

{N + n)\ AT" , 

M = ^ / - — r- 145 

N\n\ n\ ' 

The functions on fuzzy CP" have a UV cutoff given by N . Scalar fields on CP^ are 

elements in HomciVN) ■, and the integral is given by the suitably normalized trace over 

Vat. The coordinate functions Xa for a = 1, ^2 + 2?t, on fuzzy CP" are given by suitably 

rescaled generators of su{n + 1) acting on Vn- One finds [30] 

[Xa.Xb] = iANfabcX", g"^ XaXb = R"^ , (146) 

1 

dfxaXb = (n-l)(^— + -)A^x,. (147) 

n + 1 2 
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for 

R 

An = (148) 



2(n+l) ' 2 

For large A^, this reduces to the defining relations of CP" C M"^+^"'. On the other hand, 
scaling the radius as 

i?2 = N9-^ (149) 
n + 1 

near a given point (the "north pole") of CPjj gives Mg"" with U{n) invariant similarly 
as for the fuzzy sphere. We refer to [29] for further details. 

The Laplacian on fuzzy CP"' is proportional to the quadratic Casimir of su{n+l) acting 
on the functions, 

^i^) = (150) 
where Ja generates the SU{n + 1) rotations and c = It has eigenvalues 

A/,(x) = c^^^ Ux) (151) 
for fk{x) G V(fc^o,...,o,fc) according to the decomposition ()143|) . The multiplicity for given k is 

for k ^ n. To find the appropriate ao{m), consider 

[Lpn ^ _/ „ ck(k + n)/R^ + m^ ^ (n-l)\^n ^ 



2(i?/v/S)2"-i A (ky/d/R)^"-' _ 2(P/v/^)^" /-^ 



(n — 1)!2?T, ^-^ ck{k + n) / R'^ + 171"^ (n — l)\n\ Jq + rn? 

(ix^— (153) 



22"-i7r"(r2 - 1)! /n + m2 



in the large limit, where we denote the basis again with Y^"^ and used V = Vol{CP"') 
[^y^R'^^ndx = V~4 and 



This agrees with precisely with (fTT|) and generalizes the results for the fuzzy sphere. Fur- 
thermore, putting ()154|) and ()149p together gives 

/ := ^Tr(.) ^ i2ne)-Tri.) (155) 

in the above scaling limit. Therefore this rescaled fuzzy CP" is a perfect regularization of 
Mg". It has a sharp mumentum cutoff at A, and the same aQ{m) as in (j^^ can be used. 
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Appendix B: Justification of ( 1161) 



Consider first d > 3. Then 

p _ I d''x<P{xr-) Non-Planar ^ 1 /" ^ df% ,^. ^ 

and we assume that ^ m^, ^ m?. Here the integration domains have been 
rescaled to be unit balls in momentum space for all diagrams in the numerator and in 
the denominator; then the denominator can be estimated by the planar contribution which 
gives a finite contribution c (after the rescaling), which will be omitted. To proceed, let 
V oc [ki...., knd) £ K"'^ denote a unit vector in M"'^ with norm \\v\\ = 1, and consider 
generalized spherical coordinates such that d'^ki...d'^kn = dfl{v) r"''^~^dr in M"'^. Then 

d k, d k ,^.^,,ek'^ ^ / ^^^.^ / ^nd-i e ^^^^^ 



in simplified notation. The important point is that the radial integral J dr is oscillatory, 
with increasing frequency. For d > 3, the radial integral behaves as 

0(1) pOil) 

^n(d~2)-l^^A^r^vev ^ / ^n{d-2) /2-1 ^iA^uvev ^^gg) 



where = u. Clearly the (alternating) contributions increase with u. The integral is 
therefore estimated by the "last oscillation" which is of order min{l,0{l/{A'^vQv))}. We 
can exclude the region {vQv < 1/A}, whose volume goes like 0(1/ A) as A ^ oo (ignoring 
log A-corrections) . Therefore 

Rnp < 0(1 /A) (159) 

ignoring possible log-corrections^^, which establishes our claim. 

For c? = 2, we apply the same analysis to the numerator of The radial integral has 

again the form ()158|1 . but it is now dominated by small u, i.e. the "first half-oscillation". 
We therefore have to reintroduce the masses which provide an IR cutoff. To estimate this, 
we go back to the original form 

{^fd'^X^^n Non-Planar = ^ - ^.t^l ^ ^^^^'^^^ 



E fc,efc,<^ (^i)^ + m?"' {knf + rn? 

/(ln(— )). (160) 

me 



Indeed the result can only depend on the ratio — where 

me 



ml = l (161) 



15 



this is probably not the best possible estimate 
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is the NC mass scale, and in 2 dimensions it will depend only logarithmically on this ratio 
(for example, {y J d'^x(j){xY) Non-pianar = 0(ln(^)2). On the other hand, 

/ d'x<f){x)''') Planar = 0(ln(A)"), (162) 

V J m 

Hence for fixed ^ and A ^ oo, the eigenvalue distribution is again given by the Gaussian 
matrix model (jS^), however we expect the fluctuations to be larger in this case than for 
d = A. 

It is interesting to consider also the commutative scaling limit of the fuzzy sphere, 
in view of the "non-commutative anomaly" found in [24]. In that case, we can set 6 = 
1/NB? and therefore = 0{N) = 0(A). Then both planar and non-planar contributions 
have the same logarithmic behavior, and there is no well-defined eigenvalue distribution in 
that case. This was to be expected since we considered the commutative limit; however 
the "non- commutative anomaly" indicates some NC behavior in that case too, which is 
apparently too weak to induce a distinct eigenvalue distribution. The situation is different in 
4 dimensions, and a well-defined eigenvalue distribution may well exist in the commutative 
scaling limit. 
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